The Fujita phenomenon in exterior domains 
under the Robin boundary conditions 

Jean-Frangois Rault 



LMPA Joseph Liouville, ULCO, FR 2956 CNRS, 
{SJ | Universite Lille Nord de France 

50 rue F. Buisson, B.P. 699, F-62228 Calais Cedex (France) 
| | jf rault@lmpa . univ-littoral . f r| 

i> 

Abstract 

. The Fujita phenomenon for nonlinear parabolic problems dtu = Au + u p 

• in an exterior domain of K w under the Robin boundary conditions is in- 
(— i vestigated in the superlinear case. As in the case of Dirichlet boundary 

conditions (see Trans. Amer. Math. Soc 316 (1989), 595-622 and Isr. J. 
Math. 98 (1997), 141-156), it turns out that there exists a critical expo- 
nent p — 1 + 2/JV such that blow-up of positive solutions always occurs for 
subcritical exponents, whereas in the supercritical case global existence 
can occur for small non-negative initial data. 

• Key words: Nonlinear parabolic problems; Robin boundary conditions 
' Global solutions; Blow-up. 
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1 Introduction 

Let f2 be an exterior domain of , that is to say a connected open set il 
such that Q is a bounded domain when N > 2, and in dimension one, 51 is the 
■ complement of a real closed interval. We always suppose that the boundary dfl 

is of class C 2 . The outer normal unit vector field is denoted by v : 9f2 — > M. N 
and the outer normal derivative by d v . Let p be a real number with p > 1, a 
a non-negative continuous function on dQ x R + and ip a continuous function in 
Cl. Consider the following nonlinear parabolic problem 

d t u = Au + u p in £1 x (0, +oo), 

djyU -\- au == on d?l x (0, +oo), (1) 
u(-, 0) = ip in O. 

In this paper, we give a positive answer to Levine & Zhang's question pQ: the 
Fujita phenomenon, well-known in the case of f2 = (see Ref. [2]), remains 
true for the Robin boundary conditions. The case limiting a = and a = +oo 



1 



were proved by Levine & Zhang in [lj and by Bandle & Levine in [3J, respectively. 
The real number 1 + is still the critical exponent, and we prove the blowing- 
up of all positive solutions of Problem ([1]) for subcritical exponents p, whereas 
in the supercritical case, we show the existence of global positive solutions of 
Problem ([T]) for sufficiently small initial data. In the last section, we study the 
case of a general second order elliptic operator replacing the Laplacian. We also 
consider a non-linearity including a time and a space dependence. Throughout, 
we shall assume that a is non-negative 

a > on dn x R+, (2) 

and, in order to deal with classical solutions, we need some regularity on a 

a G C(dQ x R+). (3) 

To construct solutions with the truncation procedure (see Section[2]) , we suppose 

<p G C(f2), <|| ip ||oo< oo, ip > 0, lim ip(x)=0. (4) 

||x||a->oo 

In the case 17 = M. N , the boundary conditions are dropped, and the result is 
well-known by the classical paper of Fujita 2 . Thus we suppose SI ^ M. . 

2 Preliminaries 

First, we give the definition of positive solution which is understood along this 
paper. 

Definition 2.1 A positive solution of Problem (fl]) is a positive function u : 
(x, t) i->- u(x, t) of class C(U x [0,T)) nC 2A (Ti x (0,T)) ; satisfying 

d t u = Au + u p in Q x (0, +oo), 
d v u + au — on dQ x (0, +00), 
u(-, 0) = if in Q, 

where a and ip are given with (@), © and @- ^ e ti me T — T{a, ip) £ (0, +00] 
denotes the maximal existence time of the solution u. If T = +00, the solution 
is called global. 

From 3 , if T < +00, u blows up in finite time, that is to say: 

lim supu(i.t) = +00. 

Then, let us recall a standard procedure to construct solutions of Problem ([1} in 
outer domains for uniformly bounded and continuous initial data ip. For more 
details, we refer to |4] , [5] and references therein. Let (D n ) n£ ^ be a sequence 
of nested bounded domains such that 

n'cDoCfliC-C (J D n = R N . 

ties 
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Let u n be the solution of 

d t u = Au + u p in Un D n x (0, +00), 

d v u + au — on <9S1 x (0, +00), , . 

u = on dD n x (0,+oo), ^ 

u(>, 0) = y>„ in f2 n £> n , 

where ((/?n)neN denotes a sequence of functions in Cq(£1 n D n ) such that 

< ip n < tp in Q Pi D n 

and ip n —¥ ip uniformly in any compact of f2 (~l D n as n — > +00. Let z denote 
the solution of the ODE 

z = z p , 

Z(0) =|| (fi ||oo, 

with maximal existence time S — 1) ll^H 1 ' -1 " ^ ^ e corn P ar i sorl principle (see 
[B]), we have 

< u n {x,t) < u n+ i(x,t) < z(x,t) in flnD„ x [0,5]. 

Standard arguments based on a priori estimates for the heat equation imply 
u n — > u in the sense of Cf^,(fl x (0,S)) as n 4 +00, where u is a positive 
solution of Problem ([T]). Moreover, since u n vanishes on dD n for each n G N*, 
the solution u vanishes at infinity: 

lim u(i,i)=0,Vt6(0,T). 

||x||2-^O0 



3 Blow up case 

In this section, we compare the solution of Problem (fTJ) with an appropriate 
Dirichlet solution. We prove the following theorem: 

Theorem 3.1 Suppose that conditions (fj|), (® an d (0 are fullfiled. Then all 
non-trivial positive solutions of Problem ((7|) blow up infinite time for p G (1,1 + 
2/N). Moreover, if N > 3, blow up also occurs for p = 1 + 2/N. 

Proof: Ab absurdo, suppose that there exists a and a non-trivial <p satisfying 
the hypotheses above, and such that the solution u of Problem (TTJ) with these 
parameters is global. Then, consider u n the solution of the truncated Problem 
By the comparison principle from [6j, we obtain 

< u n {x, t) < u(x, t) in IT n D n for t > 0. 

Thus, u n can not blow up in finite time, and u n must be global. Next, define 
v n the solution of the following problem 

d t v n = Av n +vP in Q n D n x (0, +00), 

v„ = on dfl x (0, +00), 

v n = on dD n x (0, +00), 

«„(♦, 0) = tp n in On D„. 
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Again, the comparison principle from [5] implies < v n (x,t) < u n (x,t) in 
fi fl D n for t > 0. Then, we consider v the solution of the Dirichlet problem 

d t v = Av + v p infix (0,+oo), 
v = on dfi x (0, +oo), 

«(•, 0) = in fi, 

obtained as the limit of the v n by the truncation procedure described in Section 
[2J Thus, v < u in fix (0, +oo) and v is a global positive solution. A contradiction 
with Bandle & Levine results [3] (see [7] for the one-dimensional case). If TV > 3 
and p = 1 + 2/iV, the contradiction holds with Mochizuki & Suzuki's results [5] 
and [9]. Hence, our solution u must blow up in finite time. 



4 Global existence case 

Now, we consider supercritical exponents: 

, 2 
P>1+ N- 

We look for a global positive super-solution of Problem ([I} , we mean a function 
U satisfying 

d t u>Au + u p infix(0,+oo), 
d v u + au > on 9fi x (0, +oo), 
u(-, 0) > ip in fi. 

With this global super-solution and using the comparison principle, we construct 
the sequence (u n ) n ^ of global positive solutions of Problems ([5]). Thus, using 
the truncation procedure of Section [21 we construct a global positive solution of 
Problem ([1]). We use two different super-solutions, and we obtain two results 
on the global existence with some restrictions on the dimension N or on the 
coefficient a. First, we only suppose that the dimension 

N > 3. 

Theorem 4.1 Under hypotheses ©, (0) and @, for N > 3 and 

, 2 

Problem ([7]) admits global non-trivial positive solutions for sufficiently small 
initial data if. 

Proof: Consider ip satisfying Q and v the non-trivial positive solution v of the 
Neumann problem 

d t v = Av + v p infix(0,+oo), 
d u v = on <9fi x (0, +oo), 

v(-, 0) — (p in fi, 
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where the initial data ip is sufficiently small such that the solution v is global. 
This choice can be achieved because N > 3 and p > 1 + 2/N, see Levine & 
Zhang pp. For all a > on 9f2 x (0, +oo), we obtain 

d v v + av > on <9f2 x (0, +oo). 

Thus, w is a super-solution of Problem (frj), and we can deduce the statement of 
the theorem. 



Now, we suppose that there exists a positive constant c > such that 

a>con<90xR+. (6) 
We do not impose any condition on the dimension. 

Theorem 4.2 Let a be a coefficient satisfying (0) and @, tp an initial data 
with (0). For 

, 2 

Problem ([7]) admits global positive solutions for sufficiently small initial data ip. 
Proof: We consider the function U : f2 x [0, +oo) — > [0, +oo) defined by 

11*112 



U( X ,t)=A(t + t T»e W (-^) 



i(t + t )J' 

where [i = l/(p— 1), to > and ^4 > will be chosen below. All the calculus will 
be detailed in the proof of the general Theorem 15.31 If A > is small enough, 
we have 

d t U > AU + U p in U x (0, +oo). 
On the boundary dfl, hypothesis (HJ) gives 

d v U(x,t)+aU(x,t) > ( 2^ + to) + a ^ x ' t )) U ^ x ' t ) 

Since the boundary 90 is compact, the function (dQ 3 x 1— > — x ■ ^(x) G R is 
bounded. We choose to sufficiently big such that — x ■ v{x) / (2to) + c > 0. Then 
we obtain 

9^J7 + «[/>0onffix (0, +00). 

Finally, if we choose tp < U(-,0) in Q, the function U is a super-solution of 
Problem (TTJ). 



Remark 4.3 In the previous proof, one can note that the hypothesis (0 can be 
relaxed into 

a(x, t) > * ' "}*\ for all (x, t) G dfl x (0, +00). (7) 

Z{t + t ) 
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This condition gives us an optimal bound on a only if we know the geometry of 
the domain fl. For instance, if 

n = {\\ x \\ 2 >R}, 

we obtain 

x ■ v(x) = —R for all x £ dtt. 
Then, the equation Q is equivalent to 



a(x,t) > — for all (x,t) £ dfl x (0, +00). 

2{t + t ) 

In particular, the previous theorem holds for all non-negative a. 

In the one-dimensional case, using symmetry and translation, we can suppose 
that f2 = (—00, —1) U (1, +00). Then, without any additional hypothesis on the 
parameters of Problem ([1]), we obtain: 

Theorem 4.4 Assume the conditions (0), @ and @. For dimension N = 1 
and 

P> 3, 

Problem ((7|) admits global positive solutions for sufficiently small initial data if. 

5 Generalization 

In the manner of Bandle & Levine's results [7], we generalize our results. We 
consider the following problem 

d t u = Lu + t q || x HI vP inrix(0,+oo), 
d„u + au = on dfl x (0, +00), (8) 

u(', 0) = ip in 17, 

where q and s are two positive real numbers, p > 1 is a real number, and L 
stands for the second order elliptic operator 

N N 

L = d Xi (aij(x)d Xj ) +y^ j b i {x)d Xi . 

i,j=l i=l 

To deal with classical solutions, the coefficients are assumed to be in C 2 (f2). We 
keep the hypotheses @, (J3J) and ((U) on the parameters a and <p. In order to 
state our principal results, we shall introduce some notations. 

N 



„ » HI' 

Throughout, we assume that the matrix A = (oij)i<i,j<N is normalized, so that 
for some vq G (0, 1] 

< ^0 < P < 1 hi Ti. 
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Denote b = (pi, . . . , &jv) and let 



l(x)= ^2 [dx 3 aij{x) - 6 l (x)Jx i , 



JV 



*(x)= ^2 [p Xj aij(x) + bi(x) 



l\ 

We can state the following theorem concerning the blow-up case. 
Theorem 5.1 Assume that N > 2, 

div b(x) < in f2, 

and 

, , trace A(x) + l(x) — 
p{x) < y in Q. (9) 

Then, all non-trivial positive solutions of Problem @ 6/ow up in finite time for 

2 + 2q + s 



1 < p < 1 
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Proof: Ab absurdo, we suppose that there exists a non-trivial positive solution 
u of Problem © . As in the proof of Theorem 13.11 we deduce that there exists 
a non-trivial positive solution u of the following Dirichlet problem 

d t u = Lu + tf || x ||| v? inOx (0,+oo), 
u = on dCl x (0, +oo), 

0) = in i7. 

According to Bandle & Levine's results from [7], the solution u blows up in 
finite time under the above hypotheses. Thus, v must blow up too. 

■ 

For the one-dimensional case, Bandle & Levine weaken the hypothesis ([9]). Then, 
we obtain: 

Theorem 5.2 Assume that N = 1, 

div 6(a;) < in SI, 

and 



a + i > in Q. 



// 1 < p < 3 + 2q + s, £/ien a/Z non-trivial positive solutions of Problem Wow 
m finite time. 

Now, we consider the global existence case. 

Theorem 5.3 Assume that condition @ is satisfied, 

p(x) < 1 in O, 
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270 := inf ^ trace A + I* J > 0. 

Then, for any 

1 2 + 2q + s 
P> 1 + ^ , 

Problem @ admits global non-trivial positive solutions if the initial data if is 
sufficiently small. 



t-tt, % f / \ II x llo traceA + Z*\ rr , 



Proof: We consider the function U : £1 x [0, +00) — 5- [0, +00) denned by 

U M = A(t + t )-^, P (-^L), 

where \i = (2 + 2g + s)/(2p - 2), t > and A > will be chosen below. We 
have 

|| x III tracer 
'4(i + < ) 2 ~ 2(t + t ) 
and 

On the boundary 90, we obtain: 

d v U(x,t)+aU(x,t) = + S 

Thanks to hypothesis ([B]), and because the boundary dfl is compact, we can 
choose to sufficiently big such that 

-x ■ u{x) + c Q Qn ^ 



2*o 

Thus, d u U(x,t) + aU(x,t) > is achieved on <9£1 x (0, +00). Then, in fi, we 
have 

« TT / x t t / (, /xx II a: II? trace A + /* — 2u\ ,_. 

W *, t) - ^(x, t) = ((i - p(x)) JLJ^ + 2(t + to) " ) ^ t). 

With p < 1, we ignore the i-quadratic term, and by definition of 70, we obtain 
d t U(x,t)-LU(x,t)>(^-^)u(x,t), (10) 
with 70 — /i > 0. On the other hand, t < t + to implies 

t* || x ||; U p (x,t) < A p ~ x || x ||5 fr + to)'-**"- 1 * exp ( ~ (P 4 ~ ^ jl a "' ^(x.t). 
Using the overestimation 



(--^ 1 )%xp(^)(i + i )i >||x|||exp(- 



4(t + to) 



we obtain 



f || x HI lP{x,t) < AP- 1 ^-^)' 2 cxp^)(t + t Q )^-^-^U(x,t). (11) 

By definition of /z, we have s/2 + q — — 1) = — 1. Thus, we just have to 
choose A sufficiently small, equations (fTU)) and (ITT1) give 

d t U(x, t) - Ll/fo t) > t q \\ x \\ s 2 U p (x, t). 

Finally, if the initial data ip < U(-,0) in CI, U is a super-solution of Problem ©, 
and we can deduce the existence of a solution using the truncation procedure 
of Section [2] 
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